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Using one-dimensional spin-orbital model as a typical example of quantum spin systems with 
richer symmetries, we study the effect of an isolated impurity on its low energy dynamics in the 
gapless phase through bosonization and renormalization group methods. In the case of internal 
impurities, depending on the symmetry, the boundary fixed points can be either an open chain with 
a residual spin or (and) orbital triplet left behind, or a periodic chain. However, these two fixed 
points are indistinguishable in the sense that in both cases, the lead-correctton-to-scaling boundary 
operators (LCBO) only show Fermi-liquid like corrections to thermodynamical quantities. (Except 
the possible Curie-like contributions from the residual moments in the latter cases.) In the case of 
external (Kondo) impurities, the boundary fixed points, depending on the sign of orbital couplings, 
can be either an open chain with an isolated orbital doublet due to Kondo screening or it will fiow 
to an intermediate fixed point with the same LCBO as that of the two-channel Kondo problem. 
Comparison with the Kondo effect in one-dimensional (ID) Heisenberg spin chain and multi-band 
Hubbard models is also made. 

PACS numbers : 75.10Jm, 75.20Hr, 75.30Hx, 



I. INTRODUCTION 

In the past few years, there have been intensive 
studies on the. one-dimensionalpSpin-orbital model both 
analyticallyH^U and numericaUya^Q. Part of the reasons 
stems from the belief that the unusual magnetic prop- 
erties observed in some recently discovered qusai-one=. 
dimensional sjpin gapped materials such as Na2Ti2Sb2 0Ll 
and NaV2 05l3 can be explained by a simple two-band 
Hubbard model at quarter filling. Owing to the strong 
Coulomb repulsion, the corresponding low energy effec- 
tive Hamiltonian can then be mapped onto a quantum 
spin model: 

H = Ji Si ■ Si+i + J2 Ti ■ Ti+i 

i 

+ K {S,- S,+i){T,-T,+^) , 



J2 ~ J. The results show that when J < K/A a small 
deviation from the SU(4) point is irrelevant, hence the 
low energy properties of the model is still controlled by 
the SU(4)i fixed point mentioned above. In contrast, for 
J > K/A, the deviation results in marginally relevant 
interactions which open a gap in the spectrum, and the 
ground states are dimerized with alternating spin and 
orbital singlets. The low-lying excitations are just the 
fermions of the SO (6) Gross-Neveu model. 

Recently, the region where Ji 7^ J2 has been explored 
in Ref. and it was shown that the gapless phase can 
be extended to a large region around the original gapless 
line, (region V in Ref. || or phase B in Ref. ||). Moreover, 
the low energy physics is described by SU(2)2,s(8)SU(2)2,t 
WZNW model, and the two level-two SU(2) WZNW 
models are in general characterized by different veloci- 
(1-1) ties These conclusions are also consistent with 



where Si and Ti are spin one-half operators representing 
spin and orbital degrees of freedom at each site, respec- 
tively. For generic couplings Ji(2) and K, the model ( |l.l| ) 
has a SU(2)s ® SU(2)t symmetry. However, at the spe- 
cial couplings Ji = J2 = K /A, the symmetry group, is en- 
larged to SU(4), which is Bethe ansatz integrabled. The 
low energy effective theory at this point is known to be 
described by the SU(4)i Wess-Zumino-Novikov-Witten 
(WZNW) model, with central jiiiarge c = 3 , equivalent 
to three decoupled free bosonsH. 

Besides being as a quantum spin model for quasi-one- 
dimensional materials, the spin-orbital model can also 
appear as a low energy effective theory in other context. 
An example of this is the_spin-tube model studied re- 
cently by E.Origanc et al.EJ. In that case, the spin and 
orbital operators do not represent the real spin or orbital 
degrees of freedom, but are just mathematical objects 
used to describe the degenerate ground states obtained 
after projecting out the high energy states in the Hilbert 
space via renormalization group tran sformation. 

Earher studies of the model ( pTl| ) around the SU(4) 
point were concentrated on the Z2 symmetric case Ji = 



those obtained through numerical studieaj. 

On the other hand, the Kondo effect or more gener- 
ally speaking, the quantum impurity problem in one- 
dimensional strongly correlated electron systems is one 
of the central topics in condensed matter physics over 
the past few years. It is known that the interacting en- 
viroment developed around each local scattering center 
changes its character drastically. For example, a weak 
potential scatterer renarmalizes into an infinitely strong 
blockade to transportO, while a one-channel Kondo im- 
purity develops properties reminiscent of the two-channel 
Kondo effectEj. From theoretical point of view, these 
problems usually provide interesting realizations of non- 
Fermi liquid physics, and the advances of nanofabrication 
techniques in the last few years make many of the above 
theoretical ideas can possibly be realized in laboratories. 

In this paper, we would like to study the effect of local 
imperfection and Kondo problem of the one-dimensional 
spin-orbital model in its gapless phase. This can be viewd 
as a "stripped-down" version of quantum impurity prob- 
lems where the charge degrees of freedom have been pro- 
jected out, and may correspond to the Hubbard model 
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at special filling, which becomes insulating due to Umk- 
lapp scattering. The effects of magnetic impurities on 
ID Heisenbecg chain have already been studied by S. 
Eggert et al.tfl, and for the case of periodic chain, the 
leading-correction-to-scaling boundary operator (LCBO) 
is exactly the one that occurs in the two channel Kondo 
problem. For our case, the presence of orbital degrees 
of freedom make us expect nontrivial boundary critical 
behavior can occur through nontrivial LCBO. In fact, as 
we shall see, depending on the details, the system can 
flow to an open chain, a periodic chain or an intermedi- 
ate boundary fixed point, and nontrivial scaling behavior 
can indeed occur. The content of this paper is organized 
as follows : In section II, we discuss bosonization formu- 
las for bulk and open boundary conditions. Our results 
are presented in section III, and section IV is conclusions. 



Umklapp scatterings arising at higher order perturba- 
tion theory will result in a Mott transition at finite value 
oiU — Uc, therefore for L/ >> Uc, the charge field has 
a large gap, and only the spin-orbital part are left in the 
low energy sector. The remaining bosonized Hamiltonian 
can be further simplified by refermionization through the 
introduction of six Majorana fermions a = 1 ... 6 : 



exp{±iV'^TT'^sR{L)) , 
exp(±iV^^fRlL)) , 
exp(±iV47r$s/fl.(L)) , 



(2.4) 



where rji are Klein factors. The resulting Hamiltonian 
can then be written aaJ : 



II. BOSONIZATION FOR BULK AND 
BOUNDARY OPERATORS 



The model (|l.l|) around the SU(4) point (Ji ~ J2 ~ 
K/A) can be bosonized from the SU(4) Hubbard model 
at quarter fiUingEJ : 



U 



(2.1) 



by introducing the left and right movers for low energy 
degrees of freedom around the Fermi points {kp = Tr/Aao) 

~ Racr{x) exp{ikpx) + Laa{x) exp{—ikpx) . 



At this point, we can bosonize the above slowly varying 
fields as usual through introducing four chiral bostmic 
fields ^aaR./ L using the Abelian bosonization formulaO : 



Laiy ~ 



\/27rao 
\/27rao 



exp(iv47r$a<Tfl) , 
exp{-i\/A^<^aaL) 



(2.2) 



where the bosonic fields satisfy the commutation rela- 
tion [$a(Tfl, ^b„' l] = \5ab5„a' ' ^ud the Klein factors Kq^t 
introduced here are used to insure the anticommutation 
relations between different flavors of fermions, which sat- 
isfies the following anticommutation rule {nam k^^^'} = 
25abS„a' ■ The physical properties of the system can be 
made more transparent by changing to a new basis : 
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(2.3) 



+ {Gi + G3){ki + K2 + Kef 

+ 2G3 {ki -I- K2 + Ke){K3 + K4 + K5) 



(2.5) 



where the spin and orbital triplets are deflned as : = 
{e,e,e) and ^ {^\e,e), is defined as ^^^l- 
Gi and G2 measure the deviation from the SU(4) point, 
i.e. Ji = -J -I- Gi and J2 = ^ + G2- G3 < is a nonuni- 
versal parameter that could be extracted from the exact 
solution, and the two velocities Us,Ut are in general not 
equal to each other. 

It was shown in Ref. |^ and ^ that the Hamintonian 
( ^.5D contains several phases. Especially, there exist an 
extensive region where the system is gapless and the the 
low energy fixed point is governed by a SU(2)2,s(8)SU(2)2,t 
WZNW model : 



2tt 



Jt Jt 



(2.6) 



where ks = kt = 2, u* and are renormalized veloci- 
ties of the fixed point theory for spin and orbital sectors, 
respectively. Jg^ti^) ('^ — li2,3) are current operators 
for level two SU(2) WZNW model, whose Fourier modes 
obey the Kac-Moody algebra : 

The spin and orbital density operators have the following 
general forms: 

^ JsR + JsL + (e^^/^-^oAT, + H.c." 
+ (-l)^/""n„ 
T, ^ JtR + JtL + [e"-/^-»Mt + H.c. 

+ (-l)^/'^«nt, (2.7) 
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here J s^t are the smooth (fc ~ 0) parts of the spin( or- 
bital) density, while A^s,* and Us^t are the 2kF = 7r/2ao 
and 4fci? = 7r/ao parts. 

The current operators can be expressed in terms of 
Majorana fermions : 



i 



'tR{L) 



'2^tR(L) ^ ^tR{L) 



(2.8) 



The boson representations for 2kp components Afs,t are 

Af^ cx exp {iV7r($s + $/ + $«/)} 

- exp {iV7r(-$s + - $s/)} 
+ exp {iy/TT{<^s - */ - *s/)} 

- exp {iV7r(-$s - */ + *s/)} , 
Aft cx exp {iV7r($s + ^sf)} 

- exp - - '^sf)} 
+ exp {^^/7^(-<I>5 + - *s/)} 

- exp - + *s/)} , 
A/;+ cx exp {iy^{e, + + e,/)} 

- exp {iV7r(e^ - - ©s/)} , 
A/^^ cx exp {iv^($s + 6/ + e,j/)} 

+ exp {i + 6/ - e,/)} , (2.9) 

where 0a are the dual fields of $a, and satisfy 

[^a{x), <db{y)] = iSab<d{y - x). 

The 2kF components can be written in a more com- 
pact way by noting that the six Majorana fermions could 
be associated with six critical Ising models. Then using 
the order and disorder operators, aa and /Iq, of the Ising 
models, they can be expressed as follows : 

J\f^ cx i^J,l^J.2a■3CT4a5ae + cricr2Ai3M4M5A'6 , 
TV/ cx iala2^J-3^J.4(T5a6 + 1^^20-3 0-4^^5 Me , 

7V+ cx ((Ti^2 + i^J.l(J2) ((T3Cr4Cr5^6 + /^3Ai4M50-6) , 

7Vt+ cx (-^3a'4 + «0-3//4) im 1^2(75 fJ-6 - Cria2tJ-5<J6) ■ (2.10) 

The Akp part of spin and orbital operators, generated 
from higher harmonics of bosonization due to interac- 
tions, can be written down by noting that these operators 
should transform as vectors under S0(3)s,t and carry no 
chiralityEI : 



ris cx i^^ii A 
nt cx i^tji A 



(2.11) 



Since the fixed point is governed by a 
SU(2)2,s(8)SU(2)2,t WZNW theory, it is better to write 
the above operators in a way which makes the symme- 
try properties more transparent. This can be done by 
noting that each components of Si and Tt should trans- 
form as a vector under spin and orbital SU(2) rotations, 
respectively. This means that each components of 5'; 
and Tp^should be primary fields of the SU(2)2 WZNW 
modellia. It can then be immediately seen that Js^t are 



just the current operators of SU(2)2s,t WZNW models, 
and the 2kp components correspond to the spin 1/2 pri- 
mary fields of it. The latter can be made evident from 
Eq. ( 2.1C ) by using the equivalence between a-SU(2)2 
WZNW theory and three critical Ising modelstZl. Es- 
pecially, the spin 1/2 primary field can be expressed as 
the product of three order or disorder operators of the 
corresponding three Ising models. We then have : 



Aft 



ilsa iltO 

(2) (2) 
9s0 9ta 



^9^3(1 9to 

^9s0 9ta 



(2.12) 



where a — 1,2,3 and g^^^^^ are defined as 



(2) 



Here r" are Pauli matrices for a — 1,2,3, r° is the iden- 
tity matrix, and g^^t are spin one-half primary fields of 
SU(2)2s.t WZNW theory. The remaining primary fields 
with spin one ^[^j just correspond to the Akp compo- 
nents n-s^t- 

After completing discussions about the fixed point the- 
ory and its operator contents, we turn to the bosonized 
forms of the above operators in open boundary condition. 
The open chain boundary condition introduces the fol- 
lowing boundary ^conditions on the left- and right- mov- 
ing fermion fieldsllJ : 

Raa{0) + Laa{0) ^ , 

when transformed into boson language, it becomes 

We can then analytically continue the right-moving fields 
to left-moving fields by ^aaR{x, t) = -^-^aaL{-x, t). 
In this way, we arrive at a description of the system in 
terms of chiral fields only. 

With the above relations, we find for the boundary 
fields, we have : 



$ao-(x,i) 



'dacr{x,t) 



^$aa(0,i)--^ 



-X,t) 



*aCTL(x,<) + $a<TL(-X,i) 



2$a.L(0,i), (2.13) 



^e,,(o,i) 

or in terms of Majorana fermions : 

eR{^,t)=CL{-x.t) . (2.14) 
Substituting Eq_j2l|), Eq. (|I|) into Eq. (U|), Eq. 



9|) and Eq. ( p. 11 ), it is easy to see that all components 

of spin (orbital) operators are propotional to the current 
operators, i.e. : 

^boundary sL boundary oc JtL{0) ■ (2.15) 

This completes our discussions about the bosonization 
formulas. 
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III. BOUNDARY CRITICAL BEHAVIOR 

In this section, we shall apply the bosonization formu- 
las obtained in previous sections to discuss the possible 
boundary fixed points. Two cases are considered here : 
local defects which result in a change of local coupling 
strength compared with the bulk value and an external 
local moment coupled to the bulk system (Kondo impu- 
rity). 



A. Internal impurities 

As discussed by Eggert and AffieckEl, in the case of lo- 
cal defect, there are two important symmetries which dis- 
tinguish the possible boundary fixed points of a Heisen- 
berg chain : The site parity P5 which is refiection of the 
chain about one site, and the link parity which reflects 
the chain about one link. Another important symmetry 
of the lattice system is the translation by one site T, and 
we have the relation Pl — Ps ®T . As we shall see later, 
it is exactly the same symmetries which distinguish the 
possible boundary critical behaviors in the case of spin- 
orbital model. 

We first discuss the case where the local defect is in- 
variant under link parity Pl, i.e. altering the coupling 
strength of one link slightly. The corresponding opera- 
tors are S, -S^, T,-T,+i, and [S, ■ S,+i) x (T, -T.+i). 
Using Eq. (2T2Iand the following fusion rules for SU(2)2 



WZNW model 



Jl{z) g(uj,uj) 



z — UJ 



(3.1) 



where g and $ are spin- 1/2 and spin-1 primary fields, 
respectively. It is then easy to see that the leading con- 
tribution from Si ■ Si+i and Ti ■ Ti+i is : 

Oi — {iy {const, tr^s • tr^t -|- const, trgl ■ trgt + H.c.) . 

(3.2) 

Since both gs and gt have conformal dimensions 
(3/16,3/16) , the above operator O7 has scaling dimen- 
sion A = I and is a relevant boundary operator. Hence 
a small deviation of the coupling will be drived to strong 
coupling for either sign of SJi,2- The remaining opera- 
tor {Si ■ Si+i) X {Ti ■ Ti+i) can be extracted from fu- 
sion between two O/S and only contributes a marginal 
operator ~ tr$s + tr$t which does not affect the RG 
fiow. At this point, one important difference between 
the one-dimensional spin-orbital model and the usual an- 
tiferromagnetic Heisenberg chain should be noticed : For 
the spin-orbital model, the coupling constants between 
neighboring sites. jCap be either antiferromagnetic (AF) 
or ferromagneticaau (FM). With this in mind, we expect 
the following possible strong coupling behaviors : 



• Case I. Ji(2) > and 5Ji(2) > 0. The couplings 
will fiow toward strong AF coupling, and the local 
spin (orbital) degrees of freedom will form a singlet. 
The system becomes an open chain with two fewer 
spin (orbital) degrees of freedom. The stability of 
this strong coupling fixed point is guaranteed by 
the fact that the LCBO at this fixed point is just 
the product of two boundary spin (orbital) chiral 
current operators which have dimension two, hence 
are irrelevant. 

• Case II. Ji(2) < and 5Ji(2) < 0. The couplings 
will fiow toward strong FM couplings, and the local 
spin (orbital) degrees of freedom will form a triplet. 
The system becomes an open chain with an addi- 
tional triplet degree of freedom left. The stability 
of this strong coupling fixed point is guaranteed 
by the fact that the residual coupling between the 
triplet moment and the open chain is FM, hence is 
marginally irrelevant. 

• Case III. Ji(2) > and 5Ji(2) < or Ji(2) < 
and <5Ji(2) > 0. In this case, the coupling will flow 
to zero and leave a residual spin-orbital coupling 
K{So ■ Si){Tq ■ Ti). The fate of these degrees of 
freedom at impurity sites and 1 depends on the 
details of combinations of various possibilities. For 
example, if spin S case I, orbital S case III, the two 
sites will flrst form a spin singlet with a residual or- 
bital degrees of freedom described by ^K{Tq-Ti) 
which has a lower orbital triplet separated from a 
higher orbital singlet by a gap ~ 0{K). If in this 
case, J2 < 0, then the strong coupling fixed point is 
just an open chain with an orbital triplet. If J2 > 0, 
the triplet will be further screened by neighboring 
sites due to Kondo screening and the strong cou- 
pling fixed point is just an open chain. On the 
other hand, if spin G case II, orbital e case III, the 
two sites will form a spin triplet with a residual or- 
bital degrees of freedom described by ^K{To ■ Ti) 
which has a ground state with an orbital singlet 
separated from the higher triplet state. Then with 
the same reasoning as previous discussion, one ex- 
pect the strong coupling fixed point is either an 
open chain or an open chain with a spin triplet. In 
any case, they will still flow to either an open chain 
or an open chain with a residual ferromagnetic cou- 
pling to a spin or orbital triplet. 

The resulting possible boundary critical behavior is sum- 
marized in table |. Besides, we should mention that the 
appearance of a local triplet will not change the LCBOs, 
and only leads to an additional ground state degener- 
acy, hence an additional impurity entropy Simp = ln3. 
Of course these asymptotically decoupled local moments 
will also add a Curie- like contribution to Cimp and Ximp, 
in additional to logarithmic corrections characteristic of 
asymptotic freedom. 

We now turn to the case where the local defect re- 
spects site parity Pg, i.e. varying the coupling strength 
of two adjacent links by the same amounts. In this case. 
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the leading boundary operator arises from the sum of 
Si ■ Si+i and Ti ■ T^+i between two adjacent hnks. Due 
to the staggering factor in front of the 2kF and Akp com- 
ponents , in the continuum hmit, it is just the differential 
of Eq. dp) : 



d , 

— (const, tigs 
ax 



trgt + const. tT:g\ ■ trgt + H.c 



(3.3) 



which has scaling dimension 1 + | = |. Therefore, we 
conclude that a small deviation of coupling strength of 
two adjacent sites is irrelevant and the low energy fixed 
point is just a periodic chain. Since the open chain with 
a decoupled spin or orbital doublet is stable only for fer- 
romagnetic couplings, we arrive at the conclusion that 
when Ji.2 > 0, the open chain will be unstable and flow to 
the stable periodic chain with the impurity site included. 
However, for Ji 2 < 0, the open chain and periodic chain 
fixed points are not connected by a monotonous RG flow. 

From the above discussions, we find that for a local de- 
fect, there can be two possible boundacjy critial behaviors 
similar to that of a Heisenberg chainlU : 

• For impurities which violate the site parity Ps , the 
infrared fixed point corresponds to an open bound- 
ary condition. 

• For impurities which respect P, (hence violate Pl), 
the local defect is irrelevant if the deviation of cou- 
pling strength is not too large and at low energy, 
the chain "heals". 



TABLE I. Possible boundary critical behaviors for a local 
defect respect link parity Pl (OC represents the open chain, 
Ts and Tt represent residual spin and orbital triplets, respec- 
tively) : 





Ji > 
J2 > 


Ji > 

J2 < 


Ji < 

J2 > 


Ji < 

J2 < 


OC 








V 


OC+T, 


V 




V 


V 


OC+Tt 


V 


V 




V 


OC+Ts + Tt 











coupling is antiferromagnetic, all the above coupling con- 
stants are relevant. In fact, they satisfy the renormaliza- 
tion group equations : 



1 



dX, 



dlnL 2ttu[ 
d\B ^ 1, 
d\nL 4^ 



dX 



4kF 



1 



dlnL 



-Xf ■ X. 



4kf 



(3.5) 



Note that the most crucial difference between forward 
scattering and other couplings is that the latter breaks 
chiral symmetry and they couple orbital degrees of free- 
dom to the impurity spin. Also note that backward scat- 
tering is the most relevant one which will flow to strong 
coupling for both fcrro- and antiferromagnetic couplings. 
For the conveniences of latter boundary conformal field 
analysis, we should first transform this problem into a 
chiral form : Both the spin and orbital sectors of Hamil- 



tonian (2.6) can be mapped to a chiral theory by first 
folding the system to the positive a;— axis through defin- 



ing Jl/nix) 



EE J 



L/R 



Jl/ni^) = Jr/l{-x). Inter- 



preting the time axis as a boundary where J}^/j^(t, 0) = 
J\il{t, 0), we can then analytically continue the theory 
back to the whole a;— axis and it is described by two chi- 
ral (left-handed) currents. Furthermore, since backward 
scattering breaks the SU(2)(g)SU(2) symmetry down to 
its diagonal one, it is better to write the Hamiltonian in 
a way which makes the symmetry more transparent. r£jpr 
this purpose, we use the following coset constructionE^ : 



(3.6) 



where Q is the = 1 SUSY unitary minimal model with 
c = 1. By matching the scaling dimensions and spin 
properties, we can establish the relations between the 
product of conformal towers in these two representations. 
For our present purpose, we only need the following : 

X = (0)4 X ^2,1)] + (1)4 X [<^(2.3)] , (3.7) 



B. External impurities 

In this section, we discuss the boundary critical be- 
havior for an external spin 1/2 local moment coupled to 
the bulk system via AF exchange (Kondo) coupling. The 
Hamiltonian is decomposed as 7i = Tig + '^k where Tig 
is Eq. (p!!) and Uk is : 



'Hk — JK^irnp ■ So 



mp 
^ imp 



+ XikpTlsiO)} 



(3.4) 



Here Jk > is the Kondo coupling. Xp, Xb, X^kp are 
the forward, backward , 4fcF scattering strength, respec- 
tively, and all are proportional to Jk- Since the Kondo 



where are conformal towers of Q. The correspond- 

ing primary fields have conformal dimensions hp^q = 

{3p~2qf-l 



4g ' mI-^ ^ (^1)^ '']• ij)k is the conformal 
tower of SU(2)fe WZNW theory with spin j . With this 
new representation, the spin (orbital) current operator 
Js{0) = Jsl{0) + Jsr{0) = J\0) + J2(0) = J(0) is now 
the current operator of chiral SU(2)4 WZ^W model, and 
the Kondo interaction can be rewritten aS : 



AFJ(O) + AB*(O)0f2 



(2,3)V^(2,l) 



(3.8) 



where Xb is propotional to As, and $ is the spin one 
primary field of SU(2)4 WZNW theory. It can then be 
immediately seen that if both Xb and X4^kF vanish, .-thia. 
problem can be solved as the usual Kondo problemala-EHl 
: The impurity spin can be absorbed at special value of 
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coupling = 2u*/{2 + k)\k=A by redefining the spin 
current as the sum of that of electrons and impurity, 
J(a;) + 2'iTSimp^{x) — > J(a;). The boundary fixed point in 
this case is known to be that of the four- channel Kondo 
problem whose LCBO is the Kac-Moody descendant of 
the spin one primary field J_i • which leads to non- 
Fermi- liquid corrections to thermodynamical quantities. 
However, as we shall argue in the following, backward 
scattering will destabilize the four-channel Kondo inter- 
mediate fixed point. The system will be drived to the 
strong coupling open chain fixed point or a new inter- 
mediate fixed point depending on whether the coupling 
strength of orbital sector J2 is fcrro- or antiferromagnetic. 



In fact, from eq. (3.5), we expect that \b will scale 
to infinity first and dominate the low energy physics. To 
get a physical picture, we can do strong coupling analy- 
sis of the original lattice model : At strong coupling, the 
impurity spin forms a singlet with the spin at site zero 
and leaves an orbital doublet behaind. For the case where 
J2 < 0, the stability of this fixed point is guaranteed. We 
therefore conclude that in this case, the boundary fixed 
point is just an open chain. In other word, no LCBOs 
can come from the s=l conformal tower of the spin sec- 
tor. The leading operators are then the usual dimension 
two operators which originate from exchange interac- 
tion between edge spins and orbitals. However, in the 
case where J2 > 0, the strong coupling fixed point is 
unstable and the system will flow to some unknown in- 
termediate fixed point. To determine the fate of the RG 
flow, we consider a special point in our parameter space 
with J2 > 0, i.e. the SU(4) point where the underlying 
microscopic Hamiltonian can be considered as the SU(4) 
Hubbard model. Away from quarter filling, the Hilbert 
space can be decomposed into three parts : charge-, spin- 
and flavor conformal towers, and since the Kondo inter- 
action breaks chiral invariance, operators with nonzero 
charge, spin and flavor quantum numbers can occur as 
long as they transform as singlets under the correspond- 
ing daigonal subgroups. The possible impurity critical 
behavior can be determined by ideniihpHg LCBOs which 
satisfy the following two conditionsuEj : (i) produce a 
noninteracting limit consistent with known results, and 
(ii) respect the symmetries of the Hamiltonian. By trans- 
forming to a basis with definite parity, it is easy to see 
that only channels with positive parity coupled to the 
impurity spin, hence the noninteracting limit of SU(4) 
Hubbard model corresponds to two-channel Kondo fixed 
point. Using this as a reference point, we can borrow 
the results of Ref. O : The two leading boundary oper- 
ators are e'^^'^^'-^x (j)^ x 0/ and J^i • *^ + J?.i • *^ 
where Kp is the Luttinger liquid paremeter and 0^ is the 



charge field. 



are singlet fields under the diagonal 



SU(2) subgroup and they are equal to 0^2*1) terms of 
our previous coset language. Upon approaching quarter- 
ing, Umklapp scattering opens a charge gap for 0c and 
the remaining part of the first boundary operator has di- 
mension smaller than one, therefore should be surpressed 
by selection rules. Only the second interaction indepen- 
dent boundary operator survives in our case. Since we do 



not expect any qualitative change when the parameters 
deviate from the SU(4) point slightly as long as we are 
still in the same phase, together with the above renormal- 
ization group analysis we conclude that the system flows 
to a two channel Kondo fixed point for antiferromagnetic 
orbital couplings. 



C. Thermodynamical behavior 

In this section, we briefly discuss the corrections of 
specific heat and magnetic susceptibility induced by LC- 
BOs, i.e. the corrections of LCBO A/O(0) to the fixed 
point theory H* . The first type of LCBOs corresponding 
to open chain fixed point are the Virasoro decendants of 
identitiy operators : J^, J|, and Ji • J2. These oper- 
ators will contribute to impurity free energy defined by 
Sfunp{T,Xi) = fimpiT,Xi) - fimp{T,0) iu thc first or- 
der of A/. Since these operators have dimension two, 
we expect that A/ is proportional to l/TV, where Tfj- 
is the temperature scale at which the coupling strength 
of the relevant perturbations becomes of order one. Be- 
cause the relevant operators are of dimension 3/4, we 
have Tk oa 5J^/v for a small amount of initial change of 
coupling strength SJ . Then from dimensional consider- 
ation, we expect the correction to the specific heat Cimp 
should be proportional to T/Tk, and for the same rea- 
son, the correction to thc magnetic susceptibility Ximp 
should be proportional to 1/Tk, i.e. at T ^ 0, it pro- 
duces a T independent behavior. Note that this result 
is identical to that of the Heisenberg chain except the 
"Kondo temperature" has a different scaling relation 
with the coupling strength. 

The second type of LCBO corresponding to periodic 
chain boundary condition is Eq. (3_^) which as we shall 
see, is a Virasoro primary operator. Consequently, its 
finite-temperature expection value vanishes, and its con- 
tribution to fimp only starts from the second order of 
A/. In order to proceed the caculations, it is better to 
transform into a chiral representation as previ ous sub- 
section. To do this, we first note that Eq. (3.3) can 
be written as J_i,s • trggcrtrgt -f trgsJ-i,t ■ trgjcr, where 
J_i = Jfi,_i -I- and is equal to Ji,-i + ^2,-1 in 

terms of chiral fields. At this point, we can use the coset 



relation (3.6) to cast the above operator into a simpler 
form. In fact, noting that trgcr gf(T/^g2p 4*</'(2.3) 
and trg <7"32q ^(2,1)1 tbe leading irrelevant oper- 



ator becomes J 



s,-l 



*''<^(2,3)^(2,i) + ^ Following 
it can be shown that second-order 



methods in Ref. |2__ 
perturbation theory results in an impurity specific heat : 



1 



1 



*2A-3/4 *3/4 



„2^3-2Arp 



*2A-3/4 *3/4 



Us 



3(2A-3) 

where A = | is the dimension of -^Oj, tq is an infrared 
cutoff, and A — 3(2-1- |) — 12. Similarly, the impurity 
magnetic susceptibility can be obtained : 
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Ximp - ^2A-3/4 *3/4'^0 +0{\/T), 
Us Uf ' 

where A' = 2(2 + |)2 = 32. Note that although this 
LCBO looks nontrival, it only produces Fermi-liquid like 
behaviors because its dimension is "too high" . 

The third type of LCBO is the Kac-Moody descendent 
of spin one primary field from the spin and orbital sector 
: JL]^ ■ + J'^i ■ which appears in the case of Kondo 
impurity with antiferromagnetic orbital couplings and is 
equal to J_i • $0(3,3) in terms of coset representation, 
where J and $ are now elements oi k — A WZNW theory 
. It is exactly the same leading irrelevent operator as 
the two- channel Kondo problem except now it involves 
both spin and orbital sectors. It produces the following 
well-known fotpi of impurity specific heat and magnetic 
susceptibilityEJ : 

Ctmp ^ XjTln {Tk/T) , Xirnp ^ Xj In (Tk/T) . 



new feature in this case is that when orbital coupling 
is AF, the strong coupling fixed point will be destabi- 
lized and thermodynamical quantities at the intermedi- 
ate fixed point show interesting temperature dependence 
similar to that of the two-channel Kondo effect. 

Finally, since our theory is characterized by two differ- 
ent velocities, and in the case of periodic chain boundary 
fixed point, the contribution to Cimp and Ximp from the 
LCBO is the same as that of the dimension two Vira- 
soro descendents J^, a universal Wilson ratio can not be 
defined generally. 
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IV. CONCLUSIONS AND DISCUSSIONS 



To summarize, we have studied possible boundary crit- 
ical behaviors for the one-dimensional spin-orbital model 
in its gapless phase with a magnetic impurity. For the 
case of internal impurities, there can be either an open 
chain or periodic chain fixed point. The underlying rea- 
son for the occurrence of these two different critical be- 
haviors is similar to that of the Heisenberg spin chain 
: The leading instability for periodic chain is deter- 
mined by the spin and orbital dimerization operators 
es{x) = i^{Sj-Sj+i) andet(a;) = V{TyT-j+i) . Although 
they are allowed for impurities which violate the site par- 
ity Pg, they are prohibited for impurities respecting Ps, 
with dx£s,t{0) being the leading irrelevant operators. The 
new feature in this case is that due to the existance of 
additional orbital degrees of freedom and that the cou- 
plings between sites can be ferromagnetic, there can be 
spin or (and) orbital triplet at the impurity site together 
with the open chain at the low energy fixed point. 

For the case of Kondo impurity, we see that it can 
either flow to an open chain fixed point or an inter- 
mediate fixed point depending on the sign of orbital 
couphngs. This should be compared with the case 
of Kondo effect in Luttinger liquid or Hubbard model 
at incommensurate filling. In that case, the back- 
ward scattering which breaks chiral symmetry will re- 
sult in nontrivial leading irrelevant operators from the 
charge sector with Q_r— Ql ^ 0. The operators arising 
from coset construction obtained via diagonal embedding 
(SU(2)ij(g)SU(2)L/SU(2)diQg) can also appear together 
with the ones from the charge sector through nontrivial 
selection rules. In fact, these leading irrelevant opera- 
tors correspond to electron hopping with spin (orbital) 
flip between two ends in the open chain fixed point. In 
the present case, since the charge field is gapped, these 
processes cannot occur and only exchange interaction be- 
tween edge spins (orbitals) can appear as LCBOs. A 
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